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Abstract 

We numerically investigate quantum diffusion of an electron in a model of poly(dG)- 
poly(dC) and poly(dA)-poly(dT) DNA polymers with fluctuation of the parameters 
due to the impact of colored noise. The randomness is introduced by fluctuations 
of distance between two consecutive bases along the stacked base pairs. We demon- 
strate that in the model the decay time of the correlation can control the spread 
of the electronic wavepacket. Furthermore it is shown that in a motional narrowing 
regime the averaging over fluctuation causes ballistic propagation of the wavepacket, 
and in the adiabatic regime the electronic states are affected by localization. 

Key words: DNA, Diffusion, Correlation, Localization, Fluctuation, Motional 
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1 Introduction 

Charge transfer/transport properties in DNA attract lively interest among 
physicists, chemists and engineers p1[2|3Pf5|6f7f8] . It is now well established 
that diverse DNA structural deformations are of extreme importance during 
the charge transfer/transport process, because they help to create polarons 
which promote not only the formation of a localized electronic state, but may 
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also assist in rendering the latter mobile [9] . Remarkably, recent sophisticated 
experimental techniques, like, for example, spin-echo spectroscopy, allow to 
measure stochastic structural dynamics of monomers in polymer chains, such 
as double-stranded DNA [TU|llfl2j : There is a wealth of dynamical modes 
possessed of a quasi-continuum spectrum. In principle, each of these can influ- 
ence DNA charge transfer /transport, but, since there are more or less active 
modes [13], it is possible to take the whole manifold of DNA molecular vibra- 
tions into two parts - those which are most active, plus a "stochastic bath" 
consisting of all other ones. For the present, a number of the polaron models 
have been proposed to describe charge transfer/transport in DNA polymers, 
see, for example, Pl 1411 1 511 16fl 711 18j . On the other hand, computer simulations 
have pinpointed the crucial significance of dynamical disorder for DNA trans- 
fer/transport [20)12 Ill22|[23ll24 j. and several attempts to formulate stochastical 
models for the interplay of the former and the latter have already appeared 
in the literature, see, for example, [2^25] . 

In this communication, we shall deal with the polaron-like model by Hennig 
and coworkers as described in the works [15lll6j . where charge+breather prop- 
agation along DNA homopolynucleotide duplexes, i.e. in the both poly(dG)- 
poly(dC) and poly(dA)-poly(dT) DNA polymers has been studied. These 
works find that the coupled motion of charges and breathers connected with 
localized structural vibrations may contribute to highly efficient long-range 
conductivity. 

In our previous paper, we have investigated localization properties of electronic 
states in the adiabatic limit using a stochastic-bond-vibration approach for 
poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers within the frame- 
work of the polaron model [19]. That time, we assumed that the disorder 
is caused by DNA vibrational modes, and it influences the charge trans- 
fer/transport along DNA duplexes via electron- vibrational coupling. Here we 
present numerical results concerning the influence of stochastic changes in 
DNA hydrogen-bond stretchings and double-helical twisting angles, as well as 
the effects of finite system size, on the electron localization properties. 

Specifically, in the present paper we numerically investigate quantum diffusion 
of electron in the Hennig model of poly(dG)-poly(dC) and poly(dA)-poly(dT) 
with added fiuctuations. The latter are described by a colored noise associated 
with stochastic dynamics of the distances r{t) between two Watson-Crick base 
pair partners: < r{t)r{t ) >= rQexp(— |t — t |/r). These fiuctuations can be 
regarded as a stochastic process at high temperature, with phonon modes 
being randomly excited. In the model the characteristic decay time r of the 
correlation can control the spread of the electronic wavepacket. Interestingly, 
the white-noise limit r — can in effect correspond to a sort of motional 
narrowing regime, (see, for example, [27]) because we find that such a regime 
causes ballistic propagation of the wavepacket through homogeneous DNA 
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duplexes. Still, in the adiabatic limit r oo, DNA electronic states should 
be strongly affected by localization. 

The amplitude tq of the random fluctuations within the base pairs (the fluc- 
tuation of the distance between two bases in a base pair) and the correlation 
time r, are very critical parameters for the diffusive properties of wavepackets. 
The finding of ballistic behavior in the white-noise limit vs. localization in the 
adiabatic limit is interesting, since there is a number of experimental works 
[2S|29|3Uj observing ballistic conductance of DNA in water solutions, which is 
also temperature- independent [30] • Zalinge et al. have tried to explain the lat- 
ter effect, using a kind of acoustic phonon motions in DNA duplexes [30], which 
seems to be plausible [31], but not the only possible physical reason. We will 
propose an alternative explanation for the observed temperature-independent 
conductance, based upon our numerical results. 

The outline of the present paper is as follows. In the next section we introduce 
our DNA model for the investigation of its diffusive electronic properties. In 
Sect. 3, we give a brief explanation for the characteristic motion of wavepackets 
under the impact of colored noise. In Sec. 4 we present numerical results con- 
cerning the influence of the hydrogen-bond stretching fluctuations and twist 
angles on the localization properties. Furthermore we comment on the rela- 
tion between our numerical result and the experimental one. The last section 
contains our conclusions. 



2 Model and parameters 

The Hamiltonian for the electronic part in our DNA model is given by 

Hel = EnC^Cn - ^ (ClC'n+l + C„+iC^), (1) 

n n 

where C„ and are creation and annihilation operators of an electron at the 
site n. The on-site energies En are represented as 

E„ = Eo + A;r„, (2) 

where Eq is a constant and r„ denotes the structural fluctuation caused by the 
coupling with the transversal Watson-Crick H-bonding stretching vibration. 

The transfer integral Vnn+i depends on the three-dimensional distance dnn+i 
between adjacent stacked base pairs, labeled by n and n-|- 1, along each strand 
- and is expressed as follows. 
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Vnn+l = Vo(l - adnn+l) 



(3) 



The parameters k and a describe the strength of the interaction between the 
electronic and vibrational variables. The 3D displacements dnn+i bring about 
also a variation of the distances between neighboring bases along each strand. 
The first order Taylor expansion around the equilibrium positions is given by 



R 

dnn+i = -^(1 - cos9o){rn + r„+i). 

^0 



(4) 



Rq represents the equilibrium radius of the helix, 6^0 is the equilibrium double- 
helical twist angle between base pairs, and Eq the equilibrium distance between 
bases along one strand given by 



+ ARlsin\9o/2)) 



1/2 



(5) 



with a being the distance between neighboring base pairs in the direction of 
the helix axis. We adopt realistic values of the parameters obtained from the 
semi-empirical quantum-chemical calculations [32] • (See table 1.) 



parameter 



value 



a 

Rq 

00 

kAT 

aAT 
kcc 

O-GC 



0.1 [eF] 
0.1 [eF] 
3.4[A] 
10[A] 
36[°] 

0.778917[eyA-i] 
0.053835[A-^] 
-0.090325[eFA-^] 
0.383333[A-i] 



Table 1 

Basic parameters for DNA molecules. The subscripts, AT and GC, for k and a 
denote for ones of the poIy(dA)-poly(dT) and poly(dG)-poly(dC) DNA polymers, 
respectively. 

The Schrodinger equation describing the temporal evolution of the electron 
state vector > reads 



(6) 
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The explicit time-dependence of the Hamiltonian Hei{t) is given through both 
the time-dependence of the on-site and hopping terms. The equation can be 
expressed by scaled dimensionless variables as, 



^^eff-Q^=En{t)4>n " Kn+l(^)0n+l " K-lnW0n-l, (7) 

where 0„ =< n\ijj > and the effective Planck constant heff — 0.53. We re- 
defined the scaled dimensionless variables En{t) and Vnn+i{t) in Eq.(l), as 

En , Tp Vnn+1 , T/ 

Vb Vo ^nn+\- 

In addition to the DNA homopolymer duplexes, we also investigate the mixed 
sequence consisting of two types of the Watson-Crick pairs. Then, as a zero- 
order approximation, the electron-phonon coupling parameters for the mixed 
GC/AT stacks are taken here to be equal to the values obtained for poly(dG)- 
poly(dC) and poly(dA)-poly(dT) DNA polymers. 

We used mainly 4th order Runge-Kutta-Gill method in the numerical sim- 
ulation for the time evolution with time step 5t — 0.01. In some cases we 
confirmed the accuracy of the so obtained results by complete accord with the 
results gained with the help of a 6th order symplectic integrator that is higher 
order unitary integrations. 



3 Fluctuation of the r„(i) and motion of wavepackets 



The explicit time-dependence of the Hamiltonian comes from the fluctuation of 
the variable r„(t) at each site n. To mimic these fluctuations we use a Gaussian- 
Markovian process with standard deviation r„o (amplitude) and correlation 
time T characterized by the following covariance 



C{t -t) = < rn{t)rm{t) >= 8nmrl^e-\'-'^'\ (8) 

where < r„(t) >= 0. This process is called the Ornstein-Uhlenbeck process. 
There is no spatial correlation, viz, the stochastic fluctuations at different sites 
are independent of each other. The white-noise limit corresponds to r ^ 0. 
On the other hand, the adiabatic limit (r oo) corresponds static Anderson 
model. A numerical way for the generation of the colored noise is given in ap- 
pendix A. The fluctuations are characterized by two parameters, i.e. amplitude 
r„o and correlation time r. Another relevant quantity is. 
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oo 

Do = J C{t)dt = rl^T. 





(9) 



The quantity Dq expresses total strength for the random fluctuation and is 
related to self-diffusion coefficient [33|34|35j . The hmit r — > and r„o oo 
with keeping Dq = const, yields Gaussian white-noise characterized by < 
i^n(t)rnit ) >= ^DoSit — t). For |t — t I > r, we may treat r„(t) and r„(t ) as 
statistically independent quantities. This limiting case becomes valid when the 
lattice temperature is well above the Debye temperature. The characteristic 
decay rate Tc of the correlation function is on the order of the frequency Vt of 
the slow bond vibrations due to the relation 



re = — -0.05. (10) 

The adiabatic limit holds true for r >> allowing for the application of the 
Born-Oppenheimer approximation. 

In general, in our numerical simulation, the order of mean values < En{t) > 
,< Vrm+iit) > and the amplitude of the fluctuations (standard deviation) 
AEn{t), AVnn+lit) of the on-site energy and the transfer energy are estimated 
as follows: < En{t) >~ 0(1), AK(t) ~ fcr„o, < Kn+i(t) >~ 0(1), AV;n+i(t) ~ 
2arnoEo, where Fq = ^(1 — cos 6^0). As the fluctuating amplitude is concerned 
it can be naturally incorporated by assuming the temperature dependence as. 



rlo^ksT, (11) 

Particularly, at room temperature, i.e. for fc^T = 0.026[ey], the value of the 
scaled dimensionless thermal energy is given by ksT /Vq = 0.26. 

The motion of the electron in the fluctuating medium is crucially influenced 
by the values of the parameters r„o and r. Another important parameter is the 
band width B of the electron system which is given by the static case without 
fluctuation r„o = 0. In the scaled Hennig model for DNA the band width is 
-B ~ 2, although the exact value depends on A-T or G-C or mixed models. 
In order to infer on the onsequences of the noise for the elelctronmotion it 
illustratitive to express the parameters of the noise in units of r„o- Thus, with 
regard to its influence on the motion the noise is quantified by the effective 
parameters -B/r^o and l/(rr„o)- While B and 1/r have the effect of narrowing 
the absorption linewidth, r„o has the oppsite effect, namely broadening of 
the linewidth. The motional narrowing of the resonant absorption linewidth 
becomes important when the fluctuation rate is larger than the amplitude, 
i.e. > Tno- It then follows that time- dependent perturbation theory is 
true for B/vno » 1 or l/(rr„o) >> 1) and the adiabatic approximation is 
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applicable in the region tB » 1. In the following sections we show some 
typical wavepacket dynamics in the different parameter regimes. 

Note that the baUistic propagation of a wave packet occurs when r < 6t{— 
0.01), i.e. the fluctuation is very rapid and the system is reduced to a regular 
system without any disorder. 



4 Numerical results 



In this section, we show the numerical results of the wavepacket dynamics. 
We consider quantum diffusion of an initially localized wavepacket = 
0) = Snno- Then we mainly monitor the time-dependence of the mean square 
displacement (MSD) 

N 

m{t)=Y,<'l>n\{f>'-nof\'l>n>, (12) 
n 

and the distribution function P{n,t) = | < n|'0(i) > | — |0„(i)|. We used 
N = 2^^ and Uq = N/2 through this paper. 



4-1 Constant hopping term Kin+i = !(« = 0) 

First we present our numerical results of simple cases with a constant transfer 
integral Vnn+i = 1.0 for all n. We illustrate the typical motion of a wavepacket 
in cases when only the on-site energy En{t) fluctuates according to a Ornstcin- 
Uhlcnbcck process. Although we used mixed model in this case the qualitative 
result does not depend on the type of model even if we adapt A-T model and 
G-C model. 

Figure 1(a) shows timc-dcpcndcncc of the MSD for various correlation times 
T = 10^, 1, 0.01 at r„o = l-O- In Fig. 1(b) and (c) we depict some snap-shots of 
|0„(t)| in cases with r = 0.01 and r = 10^, respectively. In the case r = 0.01, 
■m{t) shows ballistic behavior (m(i) ~ t'^) which is due to "motional narrowing" 
caused by the short-time correlation. On the other hand, in the case r = 10^ 
m{t) shows typical localization behavior within this time scale due to Anderson 
localization. In the intermediate case t = 1, m{t) shows normal diffusion with 
Gaussian shape. 
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Fig. 1. (a) m(t) at cases with r = 0.001, 1, 100000 in the mixed model with a 
constant hopping Vnn+i = Vo = 1-0. Some snapshots of |</>n(i)l at * = 2,20,50, 100 
in the case with r = 10^ (b) and r = 0.01(c). 

4-2 Fluctuating hopping term 



Next, we show the numerical results of some cases with fluctuating hopping 
terms Vnn+i(t) obeying Eq.(3). Figure 2 shows some typical time-dependence 
of the hopping term Vi2{t) for various correlation times r = 10^, 1,0.01. The 
latter correspond to adiabatic, intermediate and rapid fluctuations, respec- 
tively. 

Figure 3 displays the time- dependence of the MSD in A-T, G-C and mixed 
model with r„o = 1-0 for various correlation times r = 10^, 1, 0.01. Snap-shots 
of |0n(^)| in A-T model for various correlation time r = 10^, 1, 0.01 are shown 
in Figure 4. As a result we obtain that the essential behavior does not depend 
on existence of the fluctuations of the hopping term. The dynamical behavior 
is qualitatively similar to the cases discussed in the preceding section with 
constant Vnn+i = 1. In the relatively short-correlation case (white-noise limit 
r = 0.01), the wave packet exhibits ballistic propagation. The extent of the 
spread of the wavepacket in the A-T model is larger than that of the G-C 
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model. On the other hand, in the adiabatic hmit (r = 10^) the wave packet 
locahzes, which corresponds to Anderson locahzation. The locahzation length 
is £ ~ 20 sites. It seems that the localization length of the A-T model is slighty 
larger than that of the G-C model on this scale. 




4000 6000 

time step 



Fig. 2. The time-dependence of the hopping term Vi2{t) at n = 2 in A-T model 
with r„o = 1.0 for various correlation times r = lO^(a), r = 1.0(b), r = 0.01(c). 

Next in Fig.5 we show a typical MSD and P{n,t) for the relatively small 
fluctuation strength r„o = 0.1 in the mixed model with r = 1. As evident 
from the relation in Eq.(9), in comparison to the case with a comparatively 
large fluctuation amphtude = 1.0 (cf. Fig. 3(b)) for a smaller amplitude 
fno = 0.1 the quantum diffusion behaves ballistically (m ~ t"^) within the same 
time scale. For the A — T and the G ~ C model we found similar behavior in 
dependence on the change in the fluctuation strength r„o. 



4-3 Diffusion rate 



Here we investigate the temporal diffusion rate defined as 



D(t) = li^^ (13) 

in the almost diffusive range 0.1 < r < 10, where < ... > expresses the average 
over several samples. If the motion of the wave packet is diffusive the diffusion 
rate is supposed to attain a constant value. Figure 6 shows the diffusion rate 
as a function of 1/r for various fluctuation amplitudes Vno in A-T and G-C 
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Fig. 3. m[t) in A-T, G-C, and mixed models with r„o = 1-0 and r = lO^(a), 
r = 1.0(b), r = 0.01(c). 

models. For the estimate of the diffusion rate we used a time interval t = lO^St. 
where the error of the diffusion rate is less than ten percent. D{t) approaches 
zero for adiabatic limit (r >> 1) due to Anderson localization, while D{t) 
goes to infinity for motional narrowing case (r << 1). As indicated in Sect. 3, 
increase of r„o means increase of temperature. It follows that the diffusion 
rate in the A-T model is larger than that in G-C model at the relatively 
low temperature (r„o = 1)- However, at higher temperatures (r„o = 3,5) the 
two models exhibit virtually equal behavior in this correlation-time range. 
Furthermore, in the high-temperature regime there seem to be no pronounced 
alterations of the diffusion coefficient as a function of r~^. 

4-4 Comparison with experimental results 

As mentioned in the introduction, in experiments concerning the temperature 
effect on the single-molecule conductance of double-stranded oligo-DNA with 
homogeneous base-pairs Zalinge et al [30j found temperature-independent bal- 
listic conductance. Furthermore, it has been shown that {dG)i5 — {dC)i5 is a 
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Fig. 4. Some snap shots of |(/>„| at t = 10,20,50 in the A-T model -with r^o = 1-0 
and T = 0.01(a), r = 1(b), r = 10^(0). 

better conductor than {dA)i5 — (ciT) 15 in the conducting experiments, just 
in accordance -with the earher experimental and theoretical findings [32j. The 
work [30] explained the observed temperature independence by deactivation 
of the acoustic modes of DNA at room and higher temperatures. 

It is interesting to recall in this context an earlier -work [36], where almost 
temperature independent photoinjected electron/hole mobility has been cal- 
culated on the basis of the stochastic Haken-Strobl-Reinecker (HSR) model 
for discotic liquid crystals. Unlike in our model, the HSR Hamiltonian used 
in [36] treats all the vibrational and lattice modes as a source of fluctuations 
of the parameters of the tight-binding model, without singling out any mode 
responsible for polaron formation. As a result, the HSR mobility temperature 
independence shows up only at small fluctuations of the molecules around 
their equilibrium positions. This is not the case in our present model. 

We have to pay attention to relatively short-time behavior and/or small spread 
regime {\/m ~ 15), as we compare the experimental result with our numerical 
results. Figure 7 shows the MSD in the A-T model with r = 0.01 and r = 
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Fig. 5. (a) m{t) and (b) \(j)n\ at t = 10,20,50 in the mixed model with 
r„o = 0.1, r = 1. 

0.0001 respectively in which cases ballistic propagation is observed due to 
the motional narrowing. In the short-time regime the time dependence of 
the MSD shows relatively weak r— dependence once the motional narrowing 
has affected the diffusive behavior. Therefore, the temperature-independence 
may be caused by the motional narrowing and the finite-size effect in the 
experiment. 

Figure 8 shows the short-time behavior of the cases r = 1 and r = 0.01 
depicted on a larger time scale in Fig.3. It follows that in the short-time 
behavior the {dG)i5 — {dC)i5 case is more diffusive than its {dA)i5 — {dT)i^ 
counterpart within the range fro which the spread of the wavepacket is -y/m ~ 
15. 



5 Summary and discussion 

We have numerically investigated quantum diffusion of an electron in the 
Hennig model of poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers 
with fluctuations of the parameters caused by colored noise. In the model 
the decay time of the correlation can control the spread of the electronic 
wavepacket. It has been shown that in a motional narrowing regime the av- 
eraging over fluctuations causes ballistic propagation of the wavepacket and 
in the adiabatic regime the electronic states are strongly localized. It has 
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Fig. 6. Diffusion rate D{t) as a function of r ^ for several fluctuation strengths 
^"710 = 1)3,5 at A-T model(a) and G-C model(b), respectively 
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Fig. 7. Short-time behavior of rn-(t) in the A-T model with r = 0.01,0.0001. 

been demonstrated that the motional narrowing affects the localization in the 
poly(dG)-poly(dC) and poly(dA)-poly(dT) DNA polymers. In either model 
temperature-dependence becomes virtually suppressed when the motion of 
the wave packet is characterized by ballistic propagation. 
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Fig. 8. Short-time behavior of m(t) of A-T, G-C, and mixed models with r„o = 1-0 
at r = 1(a) and r = 0.01(b). 



We have also investigated the temporal diffusion rate in the almost diffusive 
range. It has been found that the diffusion rate of the A-T model is larger than 
that of the G-C model at comparatively low temperatures. Interestingly for 
relatively high temperatures in the diffusive range of the wavepacket motion 
the difference between the two DNA systems gets smaller. 

Furthermore we commented on the relation between our numerical results 
and the experimental ones. It was shown that in the short-time behavior 
the significant difference of the spread of wavepackets does not exist between 
{dA)i5 — {dT)i5 and {dG)i5 ~ (c?C*)i5 except for the localization cases. 

In the present report we used periodic sequences, i.e. constant values of Eq, 
Vq, as static parts of the on-site and hopping terms for poly(dG)-poly(dC) 
and poly(dA)-poly(dT) DNA polymers, respectively. This also includes the 
mixed model. Then the motional narrowing for dynamical disorder makes the 
time-evolution of the wavepacket ballistic. However, it should be remarked 
that motional narrowing strongly localizes the wavepacket if we use disordered 
sequence for the static parts of En and/or Vnn+i- 
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A Generation of colored noise 



In this appendix, we give an algorithm in order to generate the Gaussian 
colored noise r„(m5i(:) at the mth time step ^7\. First, let us assume Gaussian 
random numbers Znim) with zero mean and unit variance at the time tm = 
m6t. Then the stochastic sequence with the colored correlation is obtained by 
the following recursion: 



where po = 0, Pm = exp( i|/r). Note that index n denotes the 

site n and m the time step m. In our numerical calculation, = const. {= 
exp(— 5t/r)). We used the algorithm independently for each site n. 

It should be noted, that there is another algorithm using the power spectrum 
of the stochastic process, which is applicable to a lot of types of the correlated 
sequence [36]l38]l3Mll35l[3iro^^ 




More recently, it has been shown that motional narrowing due to large thermal 
fluctuation effects the coherence of relaxation dynamics such as spin relaxation 
in semiconductors [27| or vibrational dephasing in a spin-Peierls system with 
lattice fluctuations [15]. 
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